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1. INTRODUCTION 

Morocco has launched a national energy strategy whose objective is to produce in 2030 at least 52% 
of electricity from renewable energy sources (RES) such as solar PV and wind energy and has progressively 
started to open up the production of electricity from RES and market it to consumers connected to the grid. 
However, the RES are intermittent and unpredictable sources because they are dependent on weather 
conditions. Thus, the integration of RES makes special requirements on the electricity grid in terms of 
quality, stability, and control of the injected power. 

The participation of Photovoltaics in energy generation is increased because of the negative 
environmental effects of conventional energy sources. In this context several publications on PV systems 
have been issued in recent years, most of the research is focused on the production side or on the grid 
injection side, or both to improve the overall efficiency. The first challenge is the extraction of maximum 
power from the PVs, in fact, several algorithms for Maximum power point tracking (MPPT) are proposed 
which can be classified into two types: classical algorithms and advanced algorithms. Among the most 
popular classical methods are the "perturb and observe" algorithm (P&O) and incremental conductance (IC) 
algorithm [1]-[3]. Farhat et al. [4] studied the effectiveness of the two techniques; (P&O) and (IC) and 
concluded that the (IC) offers a slight advantage in terms of response time dynamic. However, the measured 
power point oscillates around the MPP due to variations in the duty cycle of the DC/DC converter using a 
fixed step. Some researchers have tried to develop these methods by varying the step size [5]-[8]. Artificial 
intelligence techniques are introduced to optimize the dynamics of the PPM search such as a fuzzy logic 
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controller (FLC) [9], [10] and artificial neural networks (ANN) presented in [11], [12], these advanced 
methods have demonstrated better performance. The power injection is controlled by the grid-connected 
inverter. Hence, much research on the inverter control has been done. Noguchi et al. [13] suggests a method 
known as direct power control (DTC) which is based on the direct control of active (P) and reactive (Q) 
power without a current loop. Nevertheless, this technique has a low capacity to control the power P and Q 
simultaneously and the switching frequency is not constant [14]. To solve this problem many researchers 
have proposed various developed approaches [15]. Another technique based on indirect (P) and (Q) control 
by current control in the d-q frame is named voltage-oriented control (VOC) [16], [17]. VOC is considered to 
be more efficient due to energy loss and current distortion than DTC [16], [18]. 

This work presents the modeling and simulation of a 61 KW PV system whose objective is to 
improve the efficiency of extracting maximum power and the control of the power injected into the grid with 
a better quality of voltages and currents with a low total harmonic distortion (THD). The present study is 
structured as follows: section 2 will be devoted to the proposed method. The simulation results are detailed in 
section 3. Finally, section 4 concludes with a summary. 


2. METHOD 

To validate our approach, we based on the diagram in Figure 1. A fuzzy logic based MPPT is 
proposed to act on a boost to extract the maximum power. the VOC method with PI controllers is adopted in 
this paper for the grid side we have chosen a Three-level neutral point clamped voltage source inverter (3L- 
NPC-VSJ) controlled by the SVPWM strategy with an RL filter ensuring the elimination of harmonics. 


Boost 
Convert 3 L-NPC-inverter 


RL filter 


Figure 1. Principle diagram of the proposed control for the PV conversion system 


2.1. DC side modeling and control 
2.1.1. Photovoltaic generator (PVG) Modeling 

To set up the GPV model, we must first establish the electrical circuit equivalent. For several years 
great effort has been devoted to the study of mathematical models [9], [19]-[21]. For simplicity, the single- 
diode model as shown in Figure 2, is used in this paper. This model offers a good compromise in terms of 
ease and exactitude [9]. 

From Figure 2 it can be seen that the PV current Ipv is described by (1): 


løv = lpn — lo — lp (1) 


where: Ipv: PV generator current supplied (A), Ipn is the photocurrent that mostly depends on the solar 
irradiance and cell’s temperature, Jọ is the diode current (A) and J, is the current leak in a parallel 
resistor (A). 

The Vp-Ip characteristic of the PV array can be computed by the following (2) [9]: 
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where N, is the cell number parallel of the module, Ns is the cell number in series, Iş is the cell’s saturation 
current, q is the electron charge (1.6x10"!° C), K is the Boltzmann constant (1.38x10 J/K), Tą: is a cell 
temperature in kelvin (K), A is the ideal factor dependent on PV technology, R,: is a series resistance (Q), 
and Rp is the parallel resistor. 


Practical PV device , 


Figure 2. Equivalent model of PV cell 
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To fully understand the electrical performances of the PV array, we simulated the model using 
MATLAB/Simulink software. The characteristics of the solar panel used to validate our proposed method are 


shown in Table 1. 


Table 1. Parameters of PV module 


PV module parameter Value 
Maximum Power (W) 305.226 
Celles per module (Ncell) 96 
Open circuit voltage Voc (V) 64.2 
Short-circuit current Isc (A) 5.96 
Voltage at maximum power point Vmp (V) 54.7 
Current at maximum power point Imp (A) 5.58 
Temperature coefficient of Voc (%/deg.C) -0.2727 
Temperature coefficient of Isc (%/deg.C) 0.061745 


To get high power levels, the PV array is composed of strings of PV modules collected in parallel. 
Each string consists of modules connected in series. in our case, we used a PV array composed of 40 parallel 
strings and 5 series-connected modules per string. To understand the performance of the photovoltaic array 
we have plotted with MATLAB/Simulink as shown in Figure 3 the I-V and P-V characteristics under 
different irradiations values at a set temperature of 25 °C. From Figures 3(a) and (b), it can be seen that the 
power characteristic of the PV generator has a maximum power point MPP and the position of the MPP is 
affected by the climatic conditions (solar irradiation and temperature) of the cells. For irradiation of G = 
1kW/m? and T=25 °C, we see that the maximum power point corresponds to Pupp = 61 kW with Vinpp = 


273 V and Impp 


= 223 A. 
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Figure 3. Influence of irradiation G at temperature T=25 °C on I-V, (a) P-V and (b) characteristics 


Indonesian J Elec Eng & Comp Sci, Vol. 28, No. 1, October 2022: 76-87 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 o 79 
2.1.2. MPPT control algorithm 

In the published literature, we can find different types of algorithms allowing the search for the 
maximum power point [22]. The MPPT control must have a high level of simplicity, low power 
consumption, and low cost. The methods developed until now can be divided into two classes: artificial 
intelligence methods and conventional methods. In our study we will focus on two control algorithms: the 


classical incremental conductance (IC) [23]-[25] and proposed fuzzy logic (FL). 


a. Proposed fuzzy logic controller (FLC) 

As can be seen from Figure 4, the proposed FLC has two inputs, the error (E) and the change of 
error (CE) with the output AD which is used to vary the duty cycle D of the boost converter. The (E), (CE), 
and (AD) can be computed at sample times k by the following equations [9], [26], [27]: 


Poy (k)—Ppy(k-1) 


Et) = Vpv(K)—Vpv(k-1) = 
CE(k) = E(k) — E(k — 1) (4) 
D(k) = D(k — 1) + AD(k) (5) 


where the P,,(k),and V,,(k) are the output power and voltage of the PVG at the time k respectively. The 
FLC’s variables (E, CE, and AD) are transformed into linguistic variables. In this study, the universe of 
discourse (UOD) is divided into five fuzzy sets very low (VL), low (L), zeros (ZE), high (H), and very high 
(VH). In Figure 5 the membership functions as can be seen from Figures 5(a)-(c) triangular membership 
functions are selected for all states, except for the extremities of each function where the trapezoidal form is 
used. The inference engine is made up of a fuzzy set of rules. The development of the fuzzy rules is based on 
the use of the Mamdani strategy. 


VH 


1 NL L ZE H VL L ZE H VH 
0 
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Figure 5. The membership functions: (a) input E, (b) input CE, and (c) output AD 


To illustrate this method (Max-Min), Table 2 shows 25 rules used in our study. After combining the 
rules, the fuzzy value must be transformed into a numerical value at the output. The most popular technique 
is the centroid method, where we look for the center of gravity (CG) of the fuzzy set defined as: 
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where u( x, ) is the membership value for a point x, in (UOD), and n represents the number of rules. 


Table 2. Proposed fuzzy rule table 


CE VL L ZE H VH 
E 

VL ZE ZE VL VL VL 
L ZE ZE L L L 
ZE L ZE ZE ZE H 
H H H H ZE ZE 
VH VH VH VH ZE ZE 


2.2. AC side modelling and control 

In this section we focus on the modelling and control of a grid-connected PVG system, whose main 
objective is to study the behavior of the PV cells and to improve their performance by ensuring the extraction 
of the maximum power and the control of the inverter. The energy produced by the PV modules must be 
converted to be usable. This conversion is done through an inverter that converts DC to AC current 
synchronized with the grid, with the correct amplitude, frequency, a power factor equal 0 and minimum THD 
in conformity with the norms in order to guarantee the power quality at the point of common coupling (PCC). 


2.2.1. 3L-NPC-VSI with SVPWM 

In this paper a (3L-NPC-VSJ) shown in Figure 6 is suggested, this type of inverter offers many 
benefits compared to the classical two-level inverter 2L-VSI such as reduced voltage fluctuations (dv/dt), a 
minimizing of the THD, in addition, the semiconductors power components block a reverse voltage equal to 
only Vpc/2, and gives a compromise between the efficiency and the simplicity of control [28]-[30]. 

As can be seen from Figure 6 each arm is composed of four switches connected in series and two 
clamped diodes and each arm of the inverter can take three states. Table 3 summarizes the output voltage 
(Vio) of a 3L-NPC inverter with different switch states of phase i with i={a,b,c} and we will express these 
three different states by the notation +1,0 and -1 which correspond respectively to the voltage 
levels Vy-/2,0 and — Vac/2 .Thus the 3L-NPC-VSI has 3? = 27 possible states. Note that some switching 
states are redundant and generate the same voltage vector, so we have only 19 different voltage vectors Vp to 
Vig . This is illustrated in Figure 7. For example, the Vp vector can be generated by three different switching 
states :(1,1,1), (0,0,0), (-1, -1, -1). 


Figure 6. 3L-NPC-VSI topology connected to the grid 


Table 3. The possible switching states for each arm in 3L-NPC-VSI 
Switch states i= {a,b,c} Output voltage 


{S1iS2iS3i S3it Vio 
1100 Vicl2 ) 
0110 0 (0) 
0011 -Va/2 1) 
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Figure 7. 3-level inverter voltage vectors 


2.2.2. Vector control strategy for 3L-NPC-VSI 

The technique used in our paper is known as VOC (voltage-oriented control). This method provides 
high dynamic and static performance through internal current control loops [17]. This technique is based on 
the transformation between the stationary reference frame abc and the synchronous frame d-q using the park 
transformation. Thus, the variables are DC components in steady-state. The principle consists in measuring 
the grid voltage and detecting its angle 0 for the orientation voltage in the d q axes using the phase looked 
loop (PLL) shown in Figure 8. This loop assures the phase tracking of the direct component of the voltage, to 
eliminate the quadrature component (Eq = 0). 


Figure 8. Block diagram of the PLL used for the detection of 8 


Based on Figure 6, the inverter output voltage is calculated using the Park transformation, and 
considering that the system is balanced, the homopolar part is zero and the model is reduced to two 
components only. Considering that the d axis is in phase with the voltage grid (the Ey component and in 
phase with the voltage e,). Thus, the expressions for the d and q components of the output voltages of the 
inverter become: 


4 Rp | i d | B 
= +Lp—| P| + 7 
Val lcae æ Lg) ald ie (o 


where: Va and V} are the direct and quadratic components of the voltages at the output of the 
inverter;Ea, Eq, la, and Iq are the components of the grid voltages and currents in the d-q plane, 
respectively;w = 2mf (f = 50 Hz) corresponding to the pulsation of the three-phase grid currents; L is a 
filter inductance and R is a filter resistor. We can notice in (7), that there is a coupling between the axes of 
the park transformation, so a change on I4 provokes a change on J, and vice versa., a decoupled control for 
(7) must be applied. 


2.2.3. Indirect power control 


The control of the active and reactive powers to be injected is done through the internal loops of 
current regulation of the inverter, the relation between power and current is linear. To calculate the PI 
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correctors, the decoupled model in the park reference frame must be taken into account. To do this, we 
propose to regulate Jy and J, not by Vg and V, but by yg and y4. Thus, the following equation is obtained: 


dl 
dl 


where y4q and y4 are the output signals of the current regulators. 

It is, therefore, possible to control the I4 and J, by acting on yg and yg independently. So, it is easy 
to size the PI controller’s parameters. The expression of active and reactive powers is given in the park 
reference frame by the relation (9). 


[ = 1.5(Eqla + Eqlq) a 
Q = 15(Egla — Ealg) 


As we suppose that the d axis is synchronized with the grid voltage (Eg = Em = 220. V2 and Eq = 
0), where En is the maximum grid voltage value between phase and neutral. Thus in (10) can be simplified: 


l P = 1,5. Ea. la = Vaclac (10) 


Q = —1,5. Ea. lq 


therefore, from the relation (10), the control of the (P) and (Q) can be realized respectively by the currents I4 
and I, independently and indirectly. A PI regulator is enough for this control type. The schematic diagram of 
this control type is given in Figure 9. 


O 
vadh 


Id 


Figure 9. Diagram of the voltage-oriented control with decoupling method 


3. RESULTS AND DISCUSSION 

To validate the theoretical study and the efficiency of the proposed control strategy, a complete PV 
conversion chain structure consisting of a PVG, a boost converter, and a three-level inverter NPC connected 
to the grid is designed in the MATLAB/Simulink environment., where the parameters and data of the PV 
system module used in these simulations are given in Table 4. Two tests have been performed with different 
conditions. First, we tested the power production efficiency at different irradiation levels with a temperature 
of 25 °C and compared the two MPPT algorithms based on CI and the suggested FLC. The second test is the 
simulation of the proposed model with the same irradiation conditions using the offered fuzzy MPPT method, 
and a 3L-VSI-NPC controlled by the VOC technique with an space vector pulse width modulation 
(SVPWM) strategy. 


3.1. Testing the power production efficiency under different irradiation levels 
The simulation duration is 1.5 s and the irradiation profile adopted for our study is shown in Figure 
10 and explained by the (11). 


[Os ,0.4s]: G = 1KW/m? 
[0.45 ,0.45s]:ramp down 1 to 0.5 KW /m? 
time = [0.45 ,0.95s]:G = 0.5 KW /m? (11) 
[0.95s , 1s]: ramp up 0.5 to 1 KW /m? 
[1s ,1.5s]:G = 1 KW /m? 
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from the results presented in Figure 11, we can notice that the maximum power was extracted using the two 
algorithms MPPT CI and FLC. We see that both methods, CI and FLC, converge to the maximum power 
point MPP in steady-state. However, we observe significant oscillations around the MPP in the case of CI 
MPPT.On the other hand, the proposed FLC MPPT is stable and has negligible oscillations. This result is 
validated by Figure 12, which shows that FLC MPPT is more efficient, thus proving the effectiveness of the 
proposed technique. 


1100 + 1 r P 


Ideal Ppy 
Ppy with CI MPPT 
Ppv with FLC MPPT 


irradiation (W/m2) 
Power(KW) 


l 0.2 0.4 0.6 0.8 1 1.2 1.4 
Time (sec) S Time (sec) 


Figure 10. Irradiation profile 
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Figure 11. PV generator power response 
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Figure 12. PVG efficiency 


3.2. The control results of the grid side inverter 

As can be seen from Figure 13 we notice that the DC bus voltage response is aligned with the 
required reference (Vac*=800V) we show also a small undershoot of -10V during a ramp-down of G at 0.45s. 
We also note a neglected overshoot of 8V at 1s when a ramp-up of G is applied. Therefore, an overshoot of 
around 1.25% which is negligible, this demonstrates the performance of the PI correctors used. The 
successful response of DC voltage provides a reference current Iq* as indicated in Figure 14 which shows that 
the active current Iq is tracking their reference perfectly with a zero value of reactive current (Ig=0A). This 
means that in the steady-state, as can be seen in Figure 15, the injection of the power into the grid is done 
only with its active component (P=60.4 KW at G=1000 W/m? and P=30 KW at G=500 W/m’), while the 
reactive component of this power is zero (Q=0 KVAR). 

Figure 16 also indicates that during the steady-state, Grid voltage and current of phase A are in 
phase, which ensures a unit power factor (PF=1). Figure 17 provides the injected current, which is sinusoidal 
while being synchronized with the grid frequency (f=50 Hz). The quality of the sinusoidal current wave 
transmitted to the grid is approved by Figure 18 which presents a current THD =2.39% which conforms with 
the IEEE-519 standard (<5%). 

Figure 19 gives the output phase to phase voltage (Va) of the 3L-VSI which contains five levels of 
DC voltage (Vac, Vac/2,0, -Wac/2, -Wac). Its waveform is near the sinusoid, this voltage is filtered to reduce the 
harmonics caused by the inverter as described in Figure 20 and Figure 21. 

Figure 22 clarifies the evolution of the line-to-line voltage (Va) THD for the two topologies 2L-VSI 
and 3L-VSI-NPC with SVPWM control, it is evident that the THD in steady-state is stable at a value of about 
1.11% for the suggested method and is lower than the THD for the 2L topology which is an average of 3.3%. 
A comparative analysis between the results of the two inverter types shows that the 3L-VSI configuration 
offers higher performance while reducing the THD as much as possible. 
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Figure 13. DC bus voltage response 
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Figure 15. Active and reactive power injected into 
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Figure 17. 3-phase line currents injected into the grid 
with 3L-VSI-NPC 
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Figure 19. Line to line voltage (Va) at the 
inverter output with 3L-VSI-NPC 
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Figure 16. Grid voltage and current of phase A (pu) 
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Figure 18. current injected FFT analysis with 
proposed method 
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Table 4. Proposed system parameters 
Variable Symbol Value 
Peak output power Pmax 61 Kw 
Boost inductor Ls 5 mH 
Boost capacitor Cp 100 uF 
Boost switching frequency fs 5 Khz 
DC link Voc 800 V 
DC-Link capacitor C 1200 uF 
Inverter switching frequency fi 6 Khz 
DC controller parameters Kp, K2 20-800 
Current controller parameters Kp),Ky 0.3-20 
Filter inductor L; 2.4 mH 
Filter resistor Ry 0.001 Q 
Grid line to line voltage Vab 380 V (rms) 
Grid frequency f 50 Hz 


4. CONCLUSION 

In this paper, we have proposed a more efficient control approach for a grid-connected photovoltaic 
system under different irradiation levels. The goal of the suggested control method is to extract the maximum 
power and indirect regulation of the active and reactive power injected into the electrical network while 
regulating the DC bus voltage by ensuring that the two regulation loops (P and Q) are completely decoupled 
and independent and to inject a sinusoidal current into the network with a unit power factor and an optimal 
THD that satisfies the standard. For this purpose, an MPPT algorithm based on fuzzy logic is proposed to act 
on the boost duty cycle to obtain a higher power output. we have chosen a 3L three-phase inverter with NPC 
structure connected to the grid via an RL filter the inverter is controlled by the decoupled VOC technique 
with an SVPWM strategy. We have presented the simulation results under MATLAB/Simulink. The results 
of the different simulations prove that the suggested control model can control the power transmitted to the 
grid with good dynamics, a higher quality of the voltages, and injected currents with a minimum THD 
imposed by the standards. 
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